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Abstract 

We compute the spin-independent structure functions of the forward virtual-photon 
Compton tensor of the proton at one loop using heavy baryon chiral perturbation the- 
ory and dispersion relations. We study the relation between both approaches. We use 
these results to generalize some sum rules to virtual photon transfer momentum and 
relate them with sum rules in deep inelastic scattering. We then compute the lead- 
ing chiral term of the polarizability correction to the Lamb shift of the hydrogen and 
muonic hydrogen. We obtain — 87.05/n'^Hz and — 0.148/n^meV for the correction to 
the hydrogen and muonic hydrogen Lamb shift respectively. 



1 Introduction 



The knowledge of the forward virtual- photon Compton tensor of the proton, T'^", for energies 
of the order of the pion mass allows one to check chiral perturbation theory in the (but 
not deep) Euclidean region (Q^ = — ~ ^ 0). This is interesting by itself, since 
it allows us to test whether chiral perturbation theory computations work equally well in 
the Euclidean or Minkowski regime. Moreover, the expressions obtained for the structure 
functions permit one to generalize sum rules typically derived for real photons {q^ = 0) to 
virtual photons in the Euclidean regime (g^ ~ —"^"i < 0). These results provide us with a 
more direct connection with perturbation theory calculations, which are valid for the very 
same Euclidean sum rules but at large Q^. This idea has already been addressed in the 
case of the spin-dependent structure functions in Ref. [1]. It is the main motivation of this 
paper to obtain the spin-independent generalized sum rules. To do so we will compute the 
spin-independent structure functions at one loop in heavy baryon chiral perturbation theory 

These results are also a necessary step in the complete determination of the leading chiral- 
related polarizability corrections to the Lamb shift of the hydrogen and muonic hydroger0. 
These corrections are model independent and may yield accurate determinations of the 
proton radius provided the muonic hydrogen Lamb shift is measured with high enough 
accuracy [5J. So far only some partial results have been obtained [6J. In this paper we fill 
this gap by obtaining the complete expression for the leading correction of the polarizabity 
effect to the Lamb shift. 

We will then compute T'^" using Feynman rules and study its relation with dispersion 
relation computations in depth. We will find this computation to be a useful exercise to 
illustrate the problems and advantages of working with time-ordered correlators versus its 
non-ordered counterparts. We will clearly see that the derivation of correct dispersion rela- 
tions may turn out to be more difficult than expected, and, in particular, that the seagull 
terms are necessary in order to restore current conservation. 

The expressions for the structure functions of T'^'^ can be understood as generalized sum 
rules for ^ 0. By expanding around p = q ■ p/rnp = 0, we can relate our results with the 
sum rules in the deep Euclidean region. We will discuss the connection between the chiral 
and perturbative regime but a quantitative analysis will be relegated to future work. 

The structure of the paper is the following. In sec. 2 we define the tensors and form 
factors we use in our computation. In sec. 3 we compute W^" and T^'^ using Feynman 
rules. In sec. 4 we compute T'^" using dispersion relations and compare with the previous 
result. In sec. 5 we analyze the resulting sum rules. In sec. 6 we compute the polarizability 
corrections to the Lamb shift. In sec. 7 we give our conclusions. In the Appendix some 
contributions to T^'-' due to the A particle are computed. 

^In the same way that the spin-dependent structure functions are necessary for the hyperfine sphtting of 
the hydrogen and muonic hydrogen [3]. 
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2 Definitions 

We first define the form factors, which we will understand as pure hadronic quantities, i.e. 
without electromagnetic corrections. 

We define J^^ = J2i QiQil'^Qi where i = u,d (we will not consider the strange quark in 
this paper) and Qi is the quark charge. The form factors are then defined by the following 
equation: 



{p',s\J>'\p,s)=u{p') 



Fl(gV + ^i^2(g^)^ 



u{p) , (2.1) 



where q = p' — p and Fi, F2 are the Dirac and Pauli form factors, respectively. The states 
are normalized in the following (standard relativistic) way: 

(p', X'\p, A) = (2n)V5'(p' - p)5y, , (2.2) 

and 



u{p, s)u{p, s) = + rup) ^ , (2.3) 

where s is an arbitrary spin four vector obeying = — 1 and p ■ s — 0. 
The form factors could be (analytically) expanded as follows 

Fi{q')=Fi + ^F; + ... (2.4) 

for very low momentum. Nevertheless, we will be interested instead in their (non-analytic) 
behavior in g ~ m„. 

We also introduce the Sachs form factors: 

GE{q') = F,{q') + ^F,{q'), Guiq') = F,{q') + F,{q'). (2.5) 
We also define 

^'^^hj J'^miP^s) , (2.6) 

which has the following structure {p — q • p/rUp and we take the sign convention e°^^^ = 1): 

= + ^) W,ip, q-) + ^{p'- ^^^) - ^^^) W,{p, q-) (2.7) 

-— e^'^^'^q.s^G^ip, q') - e^^'^q,{{mpp)s, - {q ■ s)p,)G,{p, q') = W^' + W^' , 

where W^^ = W^^ and W^'^ = —W^^. Therefore, the first two (spin-independent) terms 
correspond to Wi^" and the last two (spin-dependent) terms correspond to W^^ . 
We will also need the forward virtual-photon Compton tensor 



T'^'^^iJ d^xe''i-''{p, s\T{J''{x)r{0)}\p, s) , (2.8) 
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which has the following structure {p = q ■ p/ttLj^ 
rest frame where p = q^)'- 



q although we will usually work in the 



-9^" + 2 



Slip, <l') + ^(p''- ) S,{p, q'^ 



9 \ ^ 9 / I ^ 9 

mi \ q J \ q 



e^'^P^qpS^A.ip, q') - — e^'^" qp{{m,p)s^ - (q ■ s)p^)A,{p, q'), (2.9) 



m: 



depending on four scalar functions, which we call structure functions. Similarly to W'^'^, we 
can split the tensor in the symmetric, Tg'^ = Tg^, and antisymmetric piece, T'^ = —T^. 
Moreover, it is also usual to split these tensors in two components, which we label Born 
and po/., for instance, T'^^ = T^^^^ + T^o'; • Each of these terms complies with current 
conservatioij^. The Born term is the contribution coming from the intermediate state being 
the proton (somewhat the elastic contribution) and the associated structure functions can 
be written in terms of the form factors. They read 



5f°-(p,g2) 
S^°'''ip,q') 
A'^°'\p,q') 
^f°'"'^(p,9') 



-2F^{q' 



(2mpp)2-(g2)2' 
(2mpp)2-{g2)2 ' 



-^2 \q ) + (2mpp)2-(g2)2 • 

_ 4mgpF2(g2)GM(g^) 
(2mpp)2-(g2)2 • 



(2.10) 
(2.11) 
(2.12) 
(2.13) 



From these expressions one can easily single out the point-like contributions (those due to 
a point-like particle but with non-trivial anomalous dimension, radius, ...). The remaining 
contributions, with the one-loop accuracy of our chiral computation, are encoded in the 
following expression (we split Ge,m into pieces according to its chiral counting: G 
1/m" ~ 1/A"): 



(n) 
EM 



xTr 



uu 



+ 



^(0)^(1) 



v^p+ 



qp'P+ -V p+ 



(2.14) 



qp'P^ 



where p+ = i^y^. Note that Tj^om proportional to 6{v ■ q). G^"* = 1. The expressions 
for G^^j can be found in Refs. [9l [101 lLIB- For the spin-dependent case, the only 

contribution is the term proportional to G^^G^'*, which comes from the A^™" term (this is 
the only term that contributes to the Zemach contribution to the hyperfine splitting). For 
the spin-independent case, on which we focus in this paper, we would only need G^ . 

Note that some point-like contributions are still encoded in G^^ and G'^^^ . Nevertheless, 
it should also be stressed that, in many situations, only the non-analytic behavior in the 
momentum of the structure functions is really relevant. In those cases it is not necessary to 



^The separation (definition) of the Born and pol. terms is in general ambiguous, see Refs. [71 [H! for 
an extensive discussion on this point. In our case, as far as we give an exphcit definition for T|^o^„, this 
ambiguity disappears. 



In Ref. [41 G^^ was labeled G^j , and is meant to express one loop in chiral perturbation theory. 



^(2) 
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obtain the complete expressions for the structure functions, and, in particular, counterterms 
can be avoided. This comment applies to the computation of the non-analytic behavior in 
the light quark masses and in the splitting between nucleon and the Delta mass (proportional 
to powers of l/N^ in the large N^. limit) of the Lamb shift (or hyperfine splitting) from the 
Sachs form factors. This is because only the non-analytic behavior of the structure function 
gives a non-zero contribution to the dimensional regularized one-loop computation. 
In what follows we consider the computation of , which is novel. 

3 Computation of T^^ and W^''' 

We will restrict our analysis to the SU(2) flavour case and neglect the A particle. We can 
obtain a compact expression for the polarizability contribution to W^^ at tree level (for 
simplicity we set = (1,0)): 



egA 
2A 



5(K " q") (3.15) 



_ 9!:^ , (g - 2A:)^ ) / , _ t (g - 2kr 



g° k"^ — ml + ir] } { g° k'^ — m^ -\- irj 

where K = a/p^ -fm^. A; = (0, k) = (0, q - p^), /c = g - p^. 

It is easy to check that the above expression satisfies current conservation. Prom this 
result one could try to obtain Wi and W2 directly, instead we will obtain them later from 

the imaginary part of Si. 

The expressions we obtain for 5*1 and S2 using chiral perturbation theory read 



= -(^1 m,m^ {^+'-^-{1+'-^] VT^ (3.16) 



1 


f 9A 


TT 


V2/. 




ml 







2 r o 2 / 2 

' 3 ml I . m 



q2 



,2 „0 2^, 



where 



and 



a\2 



ml 



(3.17) 



^0 



I^{mi,q\q') = / dx , ^ = (3.18) 



"TT y_ ^ rpZ 

2 — o 



q2 q2 

, , V , ( 2q^ + g^ 

— arctan ; — - + arctan 



2m^|q|; \2\q\^/ml-qi 
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gpol. 



4m2q2 + (g2) 



nip 171^, 



m 



3 /I mi ml 

- + - H 21 J 2_ 

2 + U (g0)2 



2q2 



+ 



+ 



q2 

2' 



2 - 



q- 



:i - ^) + 



, g q 



q' {q'f 



4m2q2 + (g2)2 



2 + 3^ + ^ ) Ji(m2,g0,g2^ 
q^ ' 



(3.19) 



We have obtained these results by computing the spin independent part of T*-^ = A6^^ + Bq^q^ 
from the diagrams shown in Fig. [TJ We note that only inelastic diagrams contribute to T*-^. 
T^^ has been obtained afterwards using current conservatioifl 



00 



ql 



(3.20) 



Note that this relation may miss singular terms (delta-type) in q*^ . Those are precisely the 
ones that appear in the Born term (actually current conservation holds for both the born and 
polarizability corrections independently). Polarizability terms should not have 5{q^) terms. 
From T^^ and T^, we can then reconstruct S'f°'' obtaining the results above. 





(1) 



(3) 





' P 




(2) 



I 



(4) 




/P. 

. . L 

P (Seagull) P' 

Figure 1: Diagrams contributing to T^^ . Crossed diagrams are not explicitly shown but 
calculated. 



''We have actually checked that this relation also holds by a direct computation of T*"^ using Feynman 
diagrams. 
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4 Dispersion relations 



We now try to reconstruct T'^'^ from W^" using dispersion relations. In principle, this proce- 
dure could be convenient in cliiral computations, since it reduces the number of diagrams to 
be considered, because W^'^ can be obtained with one order less. Nevertheless, in our case 
this difference is not significant. 

The structure functions of W^'^ are the ones that appear in ep inelastic scattering. We 
can relate them to the structure functions of T^*^ by using its spectral decomposition, which 
relates the imaginary term of T^*^ to W'^'^: 



Ztt 

We can use this equality to obtain Wi and W2- They read 

2 



{q\ q' 



27r 



27r2 



9a 
2U 



1 + 



m: 




qi + q^ 



V4q2(g2 



mi 



g2 + q2 _ ^4q2(g2 -^2) 



27r 2 



0{ql - ml) 



27r2 



9a 
2U 



m„ 



q 



ql 



ql 



mt 



nip nin 



q^ql 



q- 



4m2q2 + (g' 



2A2 



ql 



ql 



mt 



'mt 



q^ 



+ -W^ 2 + 3^ + ^ In 



q^ 



mt 



ql + q' + VWH 



mt 



g2 + q2 _ ^4q2(g2 



mt 



(4.21) 



mt 



(4.22) 
(4.23) 



On the other hand, by using the spectral decomposition (naively, without caring about 
convergence), we obtain the following equalitie^ 



10,th 



qo ~ql~i^' 



(4.24) 



and 



n:DR{q) = -^iT^R{q) + TUq)) 



dqo 



90, th 



qo 



ql 



(4.25) 



otherwise. It should be stressed that the equalities are derived at q^ = constant. This is quite 
different from the standard dispersion relations that are derived at = —q^ = constant. 
The latter are derived using the Cauchy theorem and the assumed analytic properties of the 
structure functions. In principle, the same results should be obtained in both approaches 
(up to, maybe, some polynomials in g°). 

Using Eq. fl4.25p . the first four (plus the crossed ones) diagrams in fig. [T]can be obtained 
using dispersion relation techniques on a diagram-by-diagram basis. If we want to obtain the 



''The expressions for T^" obtained from Feynman rules and dispersion relations may differ. Therefore, 
we label those obtained with dispersion relations with the label DR. 
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complete result using dispersion techniques the last diagram in fig. [T] is missing. This could 
be expected because, in general, the use of dispersion relations reproduces the full result up 
to, may be, some local terms, polynomials in which are usually called seagull terms. The 
last diagram can be obtained by direct computation (being a constant, it can also be fixed 
from the normalization constant of A (which is momentum independent) in the limit of real 
photons and —>■ 0, since this limit is known.). Overall we find 

= + T^aguu , (4.26) 

where 

For the "00" component we find (we do not consider S{q^) contributions, which are 
encoded in the Born term) 

T°° = T°V- (4.28) 

In this case there is not a one-to-one correspondence between the diagrams obtained from 
Feynman rules and the dispersion relation analysis. Note also that current conservation is 
lost if we only consider the result obtained from the dispersion analysis. The seagull term is 
needed for the "ij" component, whereas for the "00" component the seagull term cancels the 
no one-to-one correspondence between the Feynman diagrams and the dispersion relation 
computation. 

5 Generalized Sum Rules 

We have checked that in the limit of real photons (g^ = 0) our expressions for Sf°'" agree 
with those obtained in Ref. [9]. Obviously in the limit of g° = |q| — > 0, we also obtain the 
polarizabilities, aE and Pe, which correspond to some specific sum rules. 

We can now take the opposite limit: p — > with constant. This limit allows us a 
closer connection with the sum rules that appear in deep inelastic scattering. If we define 
the moments of deep inelastic sum rules as (we recall that x = Q'^/ {2p ■ q)) 

M^NiQ')= ['dxx''-'Wi{x,Q'), (5.29) 
Jo 

then the structure functions read (due to the symmetry of the structure functions only 
moments with even appear in our case) 

S,{p,Q')= J2 P''SI''\Q'), (5.30) 

N=0,2,i,... 

where 

HQ^) =4 r 4^W,{p,Q^) = (^-^Ym^\Q'). (5.31) 

Let us note that in the standard derivation of Eq. fl5.3ip one assumes that one can neglect 
the behavior at infinity. In our case we have an expression for Si{p, Q"^), which we can Taylor 
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expand around p = to obtain expressions for M);^(Q2) at ~ m^. Therefore, Si{p,Q^) 
by itself can be understood as a generalized sum rule, where the behavior at p ~ ~ of 
Si{p, Q^) is dictated by the integral of Wi{p', Q"^) over all values of p' weighted by the function 
l/(p'2-p2). Actually, it is argued in Ref. [1] that, in general, it is more convenient to directly 
work with Si{p^Q'^) rather than with the moments because of the existence of the inelastic 
thresholds. Either way our predictions can be tested with experiment through dispersion 
relations. Note as well that this is not necessarily equivalent to checking with experiment 
their imaginary parts at low energies, because the sum rules involve the integration over an 
arbitrary large momentum. Therefore, the naive use of Eq. fl5.3ip with the chiral expressions 
obtained in Eqs. fl4.22p and fl4.23p for Wi may yield wrong results. This is because of the 
not fast enough convergence of VFi(p, Q^) at infinity, which does not allow us to neglect its 
contribution there (nevertheless one can do so if enough derivatives are taken, but then one 
has to fix the low energy unknown constants). 

Eq. (15.291) allows us to relate in a single quantity perturbation theory and the chiral 
expansion. At large Q^, one can use perturbation theory to know their scaling with (for 
a recent state of the art see p^, for instance). At low one can use the chiral expressions 
obtained here (the combination of the elastic and inelastic terms). 

Note that in Eq. (I5.3ip both the elastic and inelastic contribution have to be added by 
definition. At large the dominant contribution comes from the inelastic term. Therefore, 
this one should match with the perturbative behavior at large Q^. On the other hand, 
our computation holds at small ~ <^ A^. In this situation the largest contribution 
may come from the elastic term. Either way one may consider sum rules for elastic or 
inelastic terms and compare them with the chiral predictions. One example could be the 
Baldin sum rule, which corresponds to the sum of the polarizabilities: aE + (^e^ and can be 
obtained from dispersion relations. The generalized Baldin sum rule would correspond (up 
to normalization) to 

(<5^) - (elastic term) . (5.32) 

For a discussion on the generalized Baldin sum rule see, for instance, Refs. [H]. In 
this paper we stop here, and relegate the quantitative analysis of the sum rules, including 
their known behavior at low (chiral) and high (perturbation theory) energies, as well as a 
comparison with experiment, to future work. 

6 Lamb shift 

An important application of our result is the evaluation of the chiral-related polarizability ef- 
fects to the Lamb shift of the hydrogen and muonic hydrogen. As we have already mentioned 
in the introduction, this result is of special importance in the case of the muonic hydrogen, 
since the measurement of the muonic hydrogen Lamb shift may yield a very precise determi- 
nation of the proton radius. At this respect polarizability effects appear as one of the main 
sources of uncertainty. Therefore, a precise (and model independent) determination of their 
value can significantly pin down the errors of the proton radius. Here we provide this value 
for the purely chiral-related effects. 

We refer to ^ (see also [^) for further details, here we only quote the main formulas 
we need for the computation. In that reference, the polarizability effects were encoded in a 
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matching coefficient named (^\fi- The expression for c^\ji can be passed to the Euchdean, 
which then reads 



3,NR 



d^kK 1 



(6.33) 



X 



{(3A;o% + k2)5i(2A;o,i^, -kl) - \^S2{iko,E, -kl)] 



where m^- is the mass of the lepton. 

Using the exphcit expressions for 5f°'- obtained in Eqs. fl3J6|) and flXTOD one could obtain 
the chiral contribution to c|']v^ due to the polarizabihty effects. In particular, it is convenient 
to rewrite c^'^vr following manner, which is amenable for a numerical analysis: 



4 2^k / 9A 



-e m. 



d'^-^k, 



f27r 



,D-1 



l + k' 



2^4 



dw 



-W 



D-5 



yj2 _^ 4-^ 1 

(1+k^)^ 



x{{2 + {l + ky)AE{w',k') + {l + k^k'w'BEiw', k') } 
where (for D = A) 



1 



dx 



1 — X 



a/1 + x'^w'^ + x{l — x)w'^k^ 



(6.34) 



(6.35) 



\ l-2x 1 /"\ (l-x)(l-2x)^ 

ax — , — - / dx- 



a/1 + x'^w'^ + x(l — x)w'^k^ "2 Jo {1 + x'^w'^ + x{l — x)w'^k'^)^ 



The functions Ae and correspond to the Euclidean version of the functions A and 
B that appear in T*-' = A6^^ + Bk^y up to a normalization factor and some momentum 
rescaling. 

Note that Eq. fl6.34p is finite. Nevertheless, one has to be careful in the computation, 
since the contributions proportional to Ae and Be are not finite by themselves and a can- 
cellation of infinities needs to take place between both terms. 

The limit of small lepton mass agrees (with logarithmic accuracy) with the result of Ref. 
[6] in Eq. (51), which could be rewritten in terms of the polarizabilities fTbj [16] . 

Following again Ref. ^ , the contribution to the potential coming from the polarizabihty 
correction reads 

5V = -^^5^^\r) , (6.36) 



and its contribution to the lamb shift reads (() stands for the average over polarizations) 

6E{n) = {E{n,l = 0) - E{n,l = 1)) = -^6io- (i^V , (6.37) 

rrip IT \ n / 

where ni^p = mimp/ {rrii^ + nip). 
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We can now put some numbers. For the case of Hydrogen we obtain 

1, 87.0488 
^^ep \a^oo = ^^■^^> 

This number should be compared with the number obtained in the logarithmic approximation 
in Ref. [6]: 

64.4841 

Hz . 



rrnpol. 



We can see that the logarithmic approximation works reasonably well. This could be ex- 
pected since the ratio me/rriT, is a small number, yet the finite piece is sizeable. In any case, 
at present, the polarizability correction is beyond the experimental accuracy of the hydrogen 
Lamb shift due to the me/m.„ suppression factor. 

For the case of the muonic hydrogen we have m^/m^r, which is not suppressed. Numeri- 
cally, we obtain 

*-rL_.^-^-V. (.30) 

In this case the logarithmic approximation does not work. On the other hand, for the muon 
case, the approximation in which the mass of the muon is equal to the pion mass in Cg'^y^ is 
a good approximation: 

T^2S'^Q 

'^^m?1a-.oo K = m^) = ^^^^V. (6.40) 

If we compare these numbers with the Zemach and vacuum polarization hadronic contri- 
butions obtained in Ref. [6], we can see that they are more or less of the same size. This 
should be contrasted with the hyperfine case [1], where the Zemach correction was the largest 
contribution (in this last case there is not vacuum polarization correction). 

Our numbers are quite similar (within the expected uncertainties) to those obtained from 
models or from dispersion relations with experimental data for the form factors [TTJ [TSj 

uniEo]. 

Finally, let us stress that our results give the leading term (in the chiral counting) to the 
Zemach and polarizability correction. Note as well that our result is parameter free, no new 
counterterm is needed. Therefore, it is model independent. The correction to our expression 
is 0{m.„/ Kqcd) suppressed. The scale of Aqcd is typically dictated by the next resonance 
that has not been integrated out, in our case the A, or, being more specific, the proton-A 
mass difference: mA — riip. This produces the largest uncertainty to our result. We expect 
to compute these corrections elsewhere. 



7 Conclusions 

We have computed the spin-independent structure functions of the forward virtual-photon 
Compton tensor of the proton at one loop using heavy baryon chiral perturbation theory 
and dispersion relations. We have studied the relation between both approaches, and used 
our results to generalize some sum rules to virtual photon transfer momentum and relate 
them with sum rules in deep inelastic scattering. We have then computed the leading 
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chiral corrections to the polarizability effect of the Lamb shift of the hydrogen and muonic 
hydrogen. We have obtained — 87.05/n^ Hz and — 0.148/n^ meV for the correction to the 
hydrogen and muonic hydrogen Lamb shift respectively. 

Extensions of this work include the computation of the A effects, as well as repeating 
the computation in SU(3) (including the strange quark). On the other hand we expect to 
perform a quantitative application of these results to the muonic hydrogen, as well as to the 
generalized sum rules. 
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A A contribution 

It is not the aim of this paper to obtain the full set of A-related contributions to the structure 
functions, which we relegate to future work, yet we can not avoid displaying a few terms 
(those which appear at tree level) that are easy to obtain. 
For the elastic A-dependent terms we obtain 



W^i(go,g') = l[^J m/(go-A)k^ 

2^ 4 / V U AM 2 



2 



W^iqo.q') = -7T hr^ m,6iqo - A)k\ (A.41) 



9 V 2mp 



16 f bi 2 m,„ A 



9 \2mpJ — — le 

16 / V / , 2 , 2n i^p^ 



A = — rUp, bf is the A— photon coupling, for extra details see P]. 

We obtain the same result for these structure functions, either if we compute them with 
dispersion relations, or with Feynman rules. 
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